Field-induced gap in the spin-i antiferromagnetic Heisenberg chain: 
A density matrix renormalization group study 
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We study the spin-i antiferromagnetic Heisenberg chain in both uniform and (perpendicular) 
staggered magnetic fields using the density-matrix renormafization-group method. This model has 
been shown earlier to describe the physics of the copper benzoate materials in magnetic field. In the 
present work, we extend the study to more general case for a systematic investigation of the field- 
induced gap and related properties of the spin-^ antiferromagnetic Heisenberg chain. In particular, 
we explore the high magnetic field regime where interesting behaviors in the field-induced gap, 
magnetization, and spin correlation functions are found. Careful examination of the low energy 
properties and magnetization reveals interesting competing effects of the staggered and uniform 
fields. The incommensurate behavior in the spin correlation functions is demonstrated and discussed 
in detail. The present work reproduces earlier results in good agreement with experimental data on 
copper benzoate and predicts new interesting field-induced features at very high magnetic field. 



PACS Numbers: 75.10.Jm, 75.40.Mg 
I. INTRODUCTION 

The quasi-one-dimensional magnetic materials have at- 
tracted considerable experimental and theoretical inter- 
est since Haldane's pioneering workld that pointed out 
the difference between the integer spin Heisenberg anti- 
ferromagnetic (HAFM) chains and the half-integer spin 
chains. By mapping the Heisenberg spin chains onto the 
0(3) nonlinear sigma model,H Haldane conjectured that 
the low-energy excitation spectrum displays a finite en- 
ergy gap for the integer spin systems while for half-integer 
spin chains it is gapless. These conjectures have been ver- 
ified by later detailed studies. 

In the linear chain HAFM family, the spin-i chain is of 
particular interest since most of its properties can be ob- 
tained exactly. These exact results serve as benchmarks 
in testing the validity of various approximation schemes. 
The low energy excitation spectrum in this system has no 
gap; the elementary excitations are spin-^ spinons; the 
ground state is quasi-long-range ordered, and the spin- 
spin correlations display power law decay. 

The effect of an applied magnetic field on the spin- 
i Heisenberg chain has also been studied to gain more 
insight into the physics of such a system. When the ex- 
ternal magnetic field is present, the Hamiltonian of the 
system is written as 



= J^(S, -Sh 



(1) 



where J is the coupling constant, Si the spin-^ operator 
on site i, Sf the z component of Si, hu = gfisH/J is 
the effective dimensionless uniform field, g the average 
effective gyromagnetic ratio and H the applied magnetic 
field. When hu = 0, the critical wave vector of the gap- 
less excitation is located at and n. The applied field 



will shift the critical wave vector of the gapless excita- 
tion away from ( for transverse spin correlations) and 
TT ( for longitudinal spin correlations) to incommensu- 
rate values while the excitation remains gapless lurtil the 
field is larger than its saturation value hu — 2.0.Q When 
hu — 2.0, the magnetization is saturated with all spins 
on the lattice sites oriented parallel to the applied uni- 
form field. Further increase in the applied field leads to 
the opening of a gap in the low energy spectrum with 
its magnitude changing almost linearly with the applied 
field, corresponding to flipping one spin to its opposite 
direction. 

The field-induced incommensurate state was first ob- 
served in the neutron scattering rneasurements on copper 
benzoate, Cu(C6D5COO)|9-3D20.|j Copper benzoate is a 
linear chain spin-i AFMq with coupling constant J ~ 
1.57 meV.El In this material, the effective spin-i Cu^+ 
ions form a linear chain structure, but the two neighbor- 
ing copper sites are not totally equivalent. Because of the 
small value of the coupling constant J, it is possible to 
study the high field (large hu) properties of Hamiltonian 
(|l|) and investigate how the induced incommensurate soft 
mode behaves with the changing magnetic field. 

In the copper benzoate experiment,0 in addition to the 
field-dependent incommensurate low energy modes, an 
unexpected non-zero energy gap induced by the mag- 
netic field was detected. The value of the gap varies 
with the magnitude and the relative orientation|-pf the 
applied magnetic fieldJd It was first suggestedH that 
the unexpected gap is caused by the inequivalence of 
the Cu sites leading to an effective staggered g ten- 
sor, which in turn gives rise to an effective staggered 
field in addition to the^japplied uniform magnetic field. 
Later, detailed analysis! shaws that an additional stag- 
gered Dzyaloshinskii-MoriyaaQ interaction term provides 
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similar contribution along with the staggered g ten- 
sor, and both have the same order in magnitude. The 
Dzyaloshinskii-Moriya interaction together with the stag- 
gered g factor can account for the observed non-zero en- 
ergy gap in copper benzoate. After making some assump- 
tions and ignoring the small exchange anisotropy, the ef- 
fective Hamiltonian to describe the ropper benzoate in 
the magnetic field can be written asiQ 

H - J^[S, • S,+i - huS^ - i-iyhsS^], (2) 

i 

where hs is the induced effective dimensionless staggered 
field, which is the key term to account for the observed 
non-zero energy gap. The magnitude of hg depends on 
the magnitude and relative direction of the applied uni- 
form field hu with respect to the sample. 

Hamiltonian (|^) has been studied using the bosoniza- 
tion approach, mapping on the Sine-Gordon model, and 
form-fact|QP--|techniques. The gap and magnetizatioe, 
behavior J3't3 the jiynamical magnetic susceptibilityjlil' L£ 
the specific health and the electron spin resonancet^ 
experiments have been analyzed using the sine-Gordon 
quantum field theory. There have also been numerical 
studies of the excitation energy and transition ampli- 
tudes in the staggered magnetic field, based on the Bethe 
Ansatz solutions.L3 These studies focus mainly on the pa- 
rameter range corresponding to tie reported copper ben- 
zoate magnetic field experiment The field- induced gap 
and related magnetization, as well as spin correlations in 
generic spin-i antiferromagnetic Heisenberg chains are of 
great interest but yet to be fully understood. Of particu- 
lar interest is the study of a wider range of parameter and 
very high magnetic field conditions beyond those probed 
by the experiment and previous theoretical work, which 
may prove to be of significant importance in understand- 
ing the entire range of field-induced phenomena and the 
underlying physics in this interesting system. 

In this paper, we report results of our numerical calcu- 
lations of the ground state and the low energy excitations 
of Hamiltonian (|^) using the density matrix renormaliza- 
tion group (DMRG) method.Ej We study the most gen- 
eral case and take hu and hs as independent variables in 
our calculations. We find that the critical (saturation) 
uniform field /i^ = 2.0 serves as an important reference 
point in understanding the obtained results. When fixed 
hu is lower than or equal to h^, the induced energy gap 
increases with hg as a power with exponent ~ |, and 
for fixed h^ larger than /ij^, the small h^ dependence is 
exponential. On the other hand, when h^ is fixed, the /i„ 
dependence of the gap displays a minimum around h^. 
When hu and hg increase simultaneously with a fixed ra- 
tio, the gap increases with the field when the ratio is small 
but develops a minimum around hu at larger ratios. The 
magnetization results are consistent with the intuitive ex- 
pectation in general. The most interesting features are 
obtained for fixed hu < h'^ and in small fixed h^ cases. 
When hu is lower than /i^ and fixed, the existence of a 



small staggered field enhances the uniform magnetization 
instead of suppressing it; similar effects are also observed 
for small fixed hg case, where the staggered magnetiza- 
tion increases until /i„ approaches hu- The uniform field 
induced incommensurate behavior is also studied, and 
the results show competing effects of the uniform and 
staggered fields, with the staggered field frustrating the 
incommensurate state. 

In the following, the numerical DMRG results will be 
presented in Section |l| and a summary given in Section 

ill. 



II. DMRG RESULTS 

The DMRG methodEl is a powerful tool for the calcula- 
tion of low lying states of quasi-one-dimensional systems 
and has been developed to calculate other properties of 
many strongly correlated systems.Ej The accuracy of the 
DMRG calculations on spin chains is generally high. This 
has also been verified in the S=l (Ref. |l^) and S=2 (Ref. 
[l^ ) Heisenberg chain calculations. In the present work, 
we employ the periodic boundary conditions (PBC) and 
use the infinite chain length algorithm of DMRG. We 
retain as many as 500 optimal states and compute up 
to chain length N=IOO in each calculation. The largest 
truncation errors are of the order of f 0~® for zero uniform 
field calculations and 10^^ for non-zero uniform field cal- 
culations. To simplify the discussion, we set the coupling 
constant J as the energy unit, hu as the effective uniform 
field and hg as the effective staggered field. 

For the Hamiltonian considered here, the effect of the 
uniform field is to induce a uniform magnetization and 
shift the critical wave vector from tt (In this paper we 
will concentrate on this case for the pitch vector, where 
a peak in the static structure factor is expected), and 
the staggered field will induce a non-zero energy gap be- 
tween the ground state and the lowest excited state. A 
non-zero staggered magnetization will also be expected 
for finite staggered fields. The nature of the ground state 
depends on the competition of the uniform and the stag- 
gered field. We will discuss the energy spectrum of the 
Hamiltonian first, and then show the ground-state mag- 
netization behavior and investigate the incommensurate 
behavior of the spin correlation functions. 

A. Energy Gap 

When the staggered field is present alone in the spin-i 
chain, that is, /iu=0 in Hamiltonian (|^), 

= ^[S,.S.+i-(-I)''/z«5f] (3) 

i 

the x-component of the total spin Sf^f. is conserved. This 
property can be used in DMRG calculations to reduce 
the dimension of the Hilbert space to be considered. 
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Thej-spin-1 case of Hamiltonian has been studied in 
detail.113 For the standard spin-1 Heisenberg chain, the 
low energy spectrum is gapful, and the lowest excited 
state is a spin triplet known as Haldane triplet. The 
presence of a staggered field will split the Haldane triplet 
into two branches, the transverse branch and the the lon- 
gitudinal branch. Both branches are gapful and the gap 
increases with the staggered field. 

For the spin-i case, the excitation spectrum for zero- 
field chain is gapless, and the presence of the staggered 
field may also open a finite gap between the ground 
state and the low-energy continuum. It is expected that 
the gapful excitations will also split into two branches, 
with the x-component of total spin Sfg^=0 (longitudinal 
branch) and 1 (transverse branch). 

The field dependence of the gap of the longitudinal and 
transverse branches for Hamiltonian (|^) is shown in Fig. 
|l| (a). It is clear that the low energy spectrum becomes 
gapful as soon as the staggered field is present. The 
magnitude of both longitudinal and transverse modes in- 
creases when the staggered field becomes larger with the 
longitudinal mode goes up faster than the transverse one. 
This behavior is exactly the same as for the spin-1 chain. 
But for the spin-1 chain in the staggered field, the in- 
crease of the longitudinal gap is about three (two) times 
faster than the transverse one for small (large) staggered 
field. For the spin-i chain case, this ratio is smaller. The 
two gaps are fitted using the equation 

A = ah'i, (4) 

where a and b are fitting parameters. The least square 
fitting gives: 

Al - flL/i^ = 2.97/^°■^^^ 

At = arh''/ = Ldlh^-^^, (5) 

where A l and A^ denote the longitudinal and transverse 
gap, respectively. The ratio of the gap increase coefficient 
of the longitudinal mode to the transverse mode ol far 
is about 1.5. The fitting curves Eq. (||) are also shown in 
Fig. |^(a). For larger staggered field, the fitting is almost 
perfect, while it is not very good for very small staggered 
field. (This deviation is not visible in Fig. |l| (a).) It 
should be noted here that the numerical error of DMRG 
results for small staggered field is much bigger than that 
for large field. 

When the uniform field in Hamiltonian (|^) is not van- 
ishing, the competition of the uniform and the staggered 
field needs to be taken into account. The presence of the 
uniform field will affect the behavior of the energy gap. 
In the non-zero uniform field case, the calculation is more 
difficult, since now even S'^^^ is not conserved any more. 
There is no good quantum number which can be used to 
reduce the relevant Hilbert space dimension. To investi- 
gate the gap behavior, the ground state and the lowest 
excited state must be calculated at the same time, and 
at least two states must be targeted in each calculation. 



To study the effect of the uniform field on the gap be- 
havior, we first calculate the staggered field dependence 
of the gap when uniform field hu is a non-zero constant. 
As described above, when the uniform field is present 
alone, there is a critical point h'^ — 2.0, where the mag- 
netization is saturated. We need to perform calculations 
at different values of hu, one in each part of the /i„ phase 
space. We choose /i„=0.5 for the small uniform field case 
and hu=3.0 for the large field case. We also calculate 
the hs dependence of the gap at the uniform field critical 
point hZ=2.0. 




1h 2 1 jj 2 3 

u u 

FIG. 1. Field dependence of the induced energy gap 
for Hamiltonian (^) with different hu and hs relations. The 
DMRG results are shown by circle, the lines denote the fitting 
curves. In (a), the filled circles denote the transverse branch, 
while the empty circles the longitudinal branch. We note that 
the "y" axis in (d) starts from A = 1 which is different with 
these starts from A = in others. 

The results for hu—0-5, 2.0 and 3.0 are shown in Fig. 
|l| (b), (c), (d), respectively. For hu—0.5 and 2.0, the 
spectrum is gapless when hs—0, and an energy gap opens 
up when the staggered field is present. The gap increases 
with the staggered field hs following the same function 
Eq. (^ as in the hu=0 case, although with different 
parameters. Least square fitting gives a=1.89, 6=0.624 
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for hu=0.5 and a=1.63, 6=0.81 for hu=2.0. Comparison 
with results in Fig. 1 (a) shows that both the coefhcient 
and the exponent in the /i„ = 0.5 case differ only slightly 
from those of the transverse branch for hu—0; for ft,„=2.0, 
the difference becomes more pronounced. The calculated 
results show that in the uniform field gapless phase < 
2.0), the induced gap is affected by the uniform field only 
when the uniform field is strong enough. For hu=3-0, 
the gap dependence is different. The energy spectrum 
is gapful even when hs=0- With the application of the 
staggered field, the gap increases. The increase of the gap 
is nearly exponential when the staggered field is not very 
strong, following A — exp{ah^) with parameters a=1.186 
and b=1.939. When the staggered field is large enough, 
the gap increase deviates from the exgpnential behavior. 

In the experiment on Cu benzoatCjcl the magnitude of 
the induced staggered field depends on the applied uni- 
form field. Roughly speaking, when the relative orienta- 
tion of the applied field is fixed, the induced staggered 
field increases linearly with the applied uniform field. To 
compare with experiment directly, we have also consid- 
ered the case when hg and /i„ increase at the same time 
with the ratio hg/hu being fixed. In Fig. |l| (e) and (f), we 
present the calculated results of the energy gap with the 
staggered field for the cases of hu — hs and hu — lOhs- 
In both cases, the field dependence of the gap for smaller 
staggered field can also be fitted with function (^). For 
hu — hg, the fitting holds at hg — hu < 2.0 with a=1.786 
and 6=0.594, while for hu = lOhg, it holds for hg < 0.1 
with a=1.755 and 6=0.613. The gap for /i„ = hg case 
increases monotonically but that for hu = lOhg exhibits 
a minimum near hg ~ 0.2 (hu 2.0), and then increases 
rapidly. The largest induced staggered field produced in 
experimeniD is about hg ^ 0.05, and it is not big enough 
to detect the gap minimum shown in Fig. |l| (f). Exper- 
iments at higher magnetic field are needed to test this 
predicted phenomenon. Here we can see again that the 
critical point /ij^ = 2.0 plays an important role in sepa- 
rating different regions where the scaling behavior of the 
field-induced gap shows qualitative difference. 

The effect of the uniform field on the spin gap induced 
by the staggered field can be studied directly by calcu- 
lating the uniform field dependence of the gap at fixed 
staggered field. From the above results, we have learned 
that the system stays in one single gapful phase when 
the staggered field is present alone, so we just choose two 
sets of hg, hg=0.05 as the small staggered field limit and 
hg=1.0 as the large staggered field limit. In Fig. (g) 
and (h) , we show the gap behavior with the change of the 
uniform field for the two cases. It is clearly seen that in 
both cases the gap decreases when the uniform field in- 
creases from zero and reaches its minimum near hu ~ 2.0, 
followed by a rapid increase with further increase in the 
uniform field. The increase of the gap after the minimum 
is approximately linear for both cases. 

Detailed calculations around the gap minimum pro- 
vide more information on the gap behavior. For hg = 0, 
the system is gapless until the uniform field reaches the 



saturation point /i^ = 2.0. When hg is not zero, the sys- 
tem has a finite energy gap, and the presence of a small 
uniform field may suppress the gap. However, for small 
uniform fields, the suppression of the gap is negligible; 
it becomes visible only when the uniform field is large 
enough. The gap minimum for non-zero hg occurs when 
the uniform field hu is slightly larger than 2.0. It occurs 
at hu ^ 2.02 for ft.s=0.05 and hu ~ 2.1 for /i^=1.0. This 
means that the minimum gap uniform field value /ij^*" 
increases slowly with hg. After the gap minimum, its 
behavior is dominated by the uniform field, and the gap 
increases almost linearly with hu- We should emphasize 
here that the minimum gap for both cases is not zero, 
and for the hg=1.0 case the decrease in magnitude is a 
little bigger than that for the hg=0.05 case. 

B. Magnetization 

The existence of a non-zero uniform field will induce a 
magnetization in the system, and the existence of a non- 
zero staggered field will induce an additional staggered 
magnetization. The staggered magnetization Mg{N) and 
the uniform magnetization M„ (TV) of the system with a 
finite chain length N are defined as 

i 

Mu{N) ^ (6) 

i 

The results for the thermodynamic hmit 
Mg = lim MgiN), 

Mu = lim Mu{N). (7) 

can be obtained by studying different chain-length sys- 
tems. The results of the magnetization and the staggered 
magnetization for the eight parameter sets {hu, hg) used 
in the previous subsection are shown in Fig. |[ 

When the uniform magnetization is absent (hu = 0), 
the induced staggered magnetization increases with the 
staggered field and approaches the saturation value 0.5 
when hg ^ oo. From the field dependence of the stag- 
gered magnetization, we can extract the staggered mag- 
netic susceptibility x*-*' = In Fig- H (a)j we can see 
clearly that the staggered magnetic susceptibility goes to 
infinity when the staggered field hg ^ instead of ap- 
proaching a constant as in the case of the spin-1 chain. 
This is because the spin-i chain is gapless for zero stag- 
gered field while it has a finite gap in the spin-1 case. 
The magnetization curves are fitted using the following 
function: 

Mg = ahl (8) 

with a=0.527 and 6=0.277. The fitting line is also shown 
in Fig. ^ (a). In the hg range shown in Fig. ^ (a), the 
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fitting is good, but it should be noted that the fitting Eq. 

will not be valid for very large hs, since it diverges 
when hs ^ oo instead of approaching the finite value 1/2 
which is the strong staggered field limit of Hamiltonian 

(i- 




FIG. 2. Magnetization curves for eight parameter sets of 
h„ and hg. The filled (empty) circles are for the staggered 
(uniform) magnetization. The solid line in (a) is the fitting 
curve described in text. For the hu~0.5 case in (b), the uni- 
form magnetization is amplified 6 times to bring it up to the 
same scale with the other panels. 

When is finite and fixed, the staggered magnetiza- 
tion increases monotonically with the staggered field and 
approaches the saturation value Ms=0.5 when hg goes to 
infinity, the same as in the hu = case. The magnetiza- 
tions for hu =0.5, 2.0, and 3.0 are shown in Fig. || (b), 
(c), and (d), respectively. Because the zero staggered 
field energy spectrum for /i„=0.5 and /i„=2.0 is gapless, 
the zero-field staggered magnetic susceptibility at hs=0 is 
still divergent. For /iu=3.0, the hg — system is gapful, 
so the zero-field staggered magnetic susceptibility has a 
finite value, x''^H0)=0.299. For /i„=2.0 and 3.0, the uni- 
form magnetization is saturated when the staggered field 
is absent; it decreases from the saturation value 0.5 when 
a staggered field is applied. When the staggered field is 



weak, the way the magnetization changes is different for 
the two cases. It decreases rapidly for /i„=2.0 but slowly 
for hu=3.0. For /i„=0.5, at zero staggered field, the uni- 
form magnetization has a finite value but is not satu- 
rated; when the staggered field increases, the uniform 
magnetization will also increase and reach a maximum 
when hs ~ 0.5 before decreasing with further increasing 
staggered field. 

In Fig. II (e) and (f), we show the field dependence of 
the magnetization for the ft-,; ~ hg and h^ — lO/ig cases, 
respectively. When the uniform and staggered fields in- 
crease simultaneously from zero, both staggered and uni- 
form magnetizations increase. For both cases, the zero- 
field staggered susceptibility is divergent and the zero- 
field uniform susceptibility is zero, we should note here 
that the zero-field uniform susceptibility is not shown 
clearly in Fig. || (e) and (f), its zero value can only be 
obtained when the results for very small field is investi- 
gated. When /i„ = hg, the effect of the staggered field 
rises rapidly and is dominating at small fields. When the 
field increases, the effect of the uniform field becomes 
more important and must be taken into account. At 
hs = ha ^ 1.2, the staggered magnetization reaches its 
maximum, and it decreases with further increasing field. 
The uniform magnetization increases monotonically with 
the increasing fields. Both staggered and uniform mag- 
netizations approach finite but nonsaturated values when 
the field goes to infinity. For /i„ — lO/ig, the effect of the 
staggered field is also dominating at small field. But 
at large field, the effect of the uniform field becomes 
dominant. The staggered magnetization decays fast to 
a small finite value 0.05 after reaching its maximum. 
The uniform magnetization increases monotonically and 
approaches a nearly saturated value 0.498 at high enough 
field. 

When the staggered field is fixed, dX hu = Q the stag- 
gered magnetization is finite while the uniform magneti- 
zation is zero. The uniform field dependence of the mag- 
netization for hs—O.Qh and hs—\.0 is shown in Fig. ||(g) 
and (h). For /is=0.05, the staggered magnetization in- 
creases with the uniform field first, reaches its maximum 
at hu ^ 1.8, and then decays to zero rapidly. For hg—l-Q, 
the staggered magnetization decreases with the uniform 
field monotonically. But the decrease is not rapid for 
small hu- In both cases, the uniform magnetization in- 
creases with the uniform field. At small uniform field, the 
uniform magnetization increases linearly with the uni- 
form field. We have obtained zero-field uniform magnetic 
susceptibility x^^HO) = 0.1185 for /is=0.05 and x^^HO) 
= 0.1283 for hs^l.Q. 

For the /i„=0.5 (/is=0.05) case, the existence of a small 
staggered (uniform) field enhances the corresponding uni- 
form (staggered) magnetization instead of suppressing 
it. This phenomenon can be explained intuitively. In 
these cases, when a small uniform (staggered) field is 
applied, the coupling between neighboring spins is weak- 
ened. While this uniform (staggered) field is not strong 
enough to destroy the effect of the stronger staggered 
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(uniform) field, it enhances the ratio between the effec- 
tive uniform (staggered) field and the effective coupling 
constant and, consequently, the induced uniform (stag- 
gered) magnetization. 



C. Correlation Function and Incommensurate 
Behavior 

We define three ground-state spin correlation functions 
for chain length L, (i) C„ parallel to the uniform magnetic 
field, (ii) Cs parallel to the staggered field, and (iii) Cy 
along the remaining {y) axis, as 



Cy{i-j)^{sysy). 



(9) 



The correlation function Cy is expected to display expo- 
nential decay because of the existence of the spin gap 
induced by the staggered field. Cs and Cu do not de- 
cay exponentially because of the effects of the non-zero 
staggered and uniform magnetization, respectively. Cu 
is also expected to show incommensurate behavior due 
to the existence of the uniform field. These correlation 
functions have the following form:! 



Cu{l) = Mt + {-iyA^ 



Vi 



■ cos(aZ -I- 6*0) 



Csii)^i-iyM^ + i-iyA 



Cy{l) 



-1)'A3- 



(10) 



where M„ and Ms are the uniform and staggered magne- 
tization, respectively, Ai, A2, A3, a, ^0 are I independent 
constants, and ^ is the correlation length. It should be 
noted here that the length ^ may not be the same along 
three different directions. 



From the spin correlation functions Eq. (p_OD, we 
can extract the corresponding correlation length. When 
hu=0, the obtained correlation length is a function of the 
staggered field. The product of the correlation length and 
the spingap gives the spin wave velocity of the system. 
In Fig. |3| we show the relation between the inverse cor- 
relation length and the transverse gap. It is seen that 
the curve goes linearly which means that the spin wave 
velocity does not change with the magnitude of the stag- 
gered field. The linear fitting of the line gives the spin 
wave velocity v — At^t = 1/0.6364 — 1.5713, in good 
agreement with the exact spin wave velocity for thfi-spin- 
i chain ^ « 1.5708. As in the spin-1 chain case,Ej this 
also serves as an independent check of the self-consistency 
of our calculations. From the figure, we can see for zero 
stagger field spin-i chain, the spectrum is gapless, A = 0, 
so the correlation length is infinite. For other cases in our 
calculations, the correlation length is difficult to obtain, 
because the chain length is short and the numerical er- 
ror is bigger for those cases. But we can conclude from 
our results that in general the spin wave velocity is not 
a constant any more, and, instead, it changes with the 
applied uniform field. 
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FIG. 4. Incom- 

mensurate behavior in = \/le'^^{Cu{l) — M^) for hu=0-5, 
hs— 0.05 and chain length L=100 system. The correlation 
length ^ « 4.05. The odd and even / are denoted by filled and 
empty circles, respectively. The solid and dashed fitting lines 
are ±0.08 cos(0.38Z - 0.407). 
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FIG. 3. The inverse correlation length for the trans- 
verse branch of the zero uniform field case vs the transverse 
gap at difi^erent staggered fields. The solid fitting line is 
= 0.6364 At. 



The cosine function in C„ comes from the incommen- 
surate behavior induced by the uniform field. Because of 
the existence of the staggered field, the net incommensu- 
rate behavior is a result of competition between and 
hs. In Fig. |, we present 0(/) = yfle'/^Cuil) - Ml) 
at chain length 100 as a function of / for /iu=0.5 and 
/is=0.05, where M„ is obtained by the magnetization 
calculation discussed above. This clearly shows the exis- 
tence of incommensurability in the system. In this case, 
the correlation length ^ « 4.05. 
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FIG. 5. Static structure factor Su{q) for /iu=0.5 and 
/is=0.05. The results for chain length from 60 to 100 are 
shown. 

The incommensurate behavior in the correlation func- 
tion leads directly to the peak shift from tt in the static 
structure factor S{q) which can be obtained from the cor- 
relation functions. For the uniform field (z) axis and the 
staggered field {x) axis, we can write: 



(11) 
(12) 



Since we use the periodic boundary conditions in our cal- 
culations, the wave vector q is well defined, q = n = 
1,...,L. 

In Fig. H, we present the static structure factor 5„ at 
hu = 0.5 and hs—Q.Ob for the even chain length from 60 
to 100. The results for different chain length systems fall 
onto the same curve. This success is due to the small 
correlation length of the system considered ~ 4.05). 
The chain lengths used here are much larger than the 
correlation length, and the finite size effect is very weak. 
Fig. 1^ shows a two peak structure symmetric about tt 
which is obtained by Eq. (jl^). For this case, we have 
obtained the critical wave vector shift 5q — \q — 'k\ ^ 
0.224 ^ 0.07 TT. The accuracy of 5q obtained from the 
peak deviation of Su{q) is good; it is mainly limited by 
the finite chain lengths used in the calculations and the 
error is estimated to be less than 1 %. 
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FIG. 6. The field dependence of the critical wave vector 
shift 5q for (a) ftu=0.5; (b) /i„ — Whs; and (c) = 0.05. 

For all parameters studied in these calculations, using 
the peak position of the static structure factor Su{q) to 
determine the existence of the incommensurate state, we 
found the critical wave vector shift in three sets of pa- 
rameters, hu ~ 0.5, hu — lOhs, and hg ~ 0.05. The 
critical wave vector shift 6q versus field in these cases is 
shown in Fig. ^. Since the largest total chain length in 
our calculations is N=100, we cannot detect the incom- 
mensurate state if Sq < = 0.027r. For the hu = 0.5 
case, the critical wave vector shift becomes smaller when 
the staggered field hg increases from zero. It is expected 
to go to zero when the staggered field is large enough. In 
other words, the existence of the staggered field will frus- 
trate the incommensurate state. For the other two cases 
hu = lOhs and kg = 0.05, the critical wave vector shifts 
away from tt with increasing uniform field but eventually 
returns to tt around /i^. 

For comparison, we examine our calculated results 
with parameters corresponding to the reported experi- 
mental case on Cu benzoateJj The highest magnetic field 
reached in the experiment is about 7T which corresponds 
to hu ~ 0.52 . The energy gap observed at this field is 
A ~ 0.4 meV, considering the coupling constant for the 
material J=1.57 meV, A/ J ~ 0.255. From our calcu- 
lated results for /i„ = 0.5, we estimate that the induced 
staggered field hg ~ 0.04. The critical wave vector shift 
5q for hu^O.5, and hg=OM is 6q - 0.274 - 0.09 tt. The 
critical wave vector shift for hu=0.52, /is=0.04 will be 
slightly larger than 0.09 tt. In the experiment, the largest 
wave shift for H=7T is about 0.12 tt. This comparison 
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shows that the calculated results are in good agreement 
with the experiment. 

III. SUMMARY 

We have carried out systematic calculations using the 
density matrix renormalization group method to study 
the behavior of the energy gap, magnetization, and spin 
correlation functions of spin-i antiferromagnetic Heisen- 
berg chain in the presence of a uniform and a perpen- 
dicular staggered magnetic field. An extensive examina- 
tion of the parameter space has revealed many interest- 
ing features beyond those reported in previous studies. 
In particular, results at very high magnetic field show 
quantitatively and even qualitatively different behaviors 
in the energy gap and magnetization, from those found 
for lower fields. For the = lOhs case, which is close to 
the real parameters in the experiments, we can see from 
Fig. 1^ and Fig. |^ that the competition of the staggered 
and uniform field is visible only when hg > 0.15, that is 
hu > 1.5. In Cu benzoate, hu = 1.5 corresponds to a 
applied field ~ 21 Tesla. Further experimental inves- 
tigation at magnetic fields higher than 21 Tesla is needed 
to test these predictions. 

The field-induced energy gap is dominated by the stag- 
gered field when the uniform field is smaller than the 
standard spin-i chain saturation field /i^=2.0. When 
the uniform field is comparable or larger than /i^, the 
effect of the uniform field becomes important and must 
be taken into account. The uniform field introduces frus- 
tration effects and creates a local minimum in the energy 
gap near h'^ in several cases. 

The magnetization results clearly reveal the competi- 
tion of the uniform and staggered fields. However, for 
some uniform (staggered) field, the existence of a small 
staggered (uniform) field enhances the uniform (stag- 
gered) magnetization instead of suppressing it. The com- 
petition of the two kinds of field can also be seen from 
the incommensurate behavior with the staggered field 
suppressing the incommensurate state and moving the 
critical wave vector closer to the zero-field valuCpSr. 

The specific heat experiment on Cu benzoatg3 shows 
that the field-induced gap scales with approximately 
2/3 power jn£.the applied magnetic field. The ana- 
lytic resultsD^EiJ yield the same power law dependence for 
Hamiltonian (|^). In our numerical results, for — 
case, the scaling power 6^ =0.678 for the longitudinal 
gap and 6t=0.63 for the transverse gap. They are in 
good agreement with the experiment and the analytic 
value. The existence of the non-zero uniform field mod- 
ifies the power law relation. When the uniform field is 
small, the modification is almost negligible. Further high 
field experiments are needed to examine the predicted ef- 
fect of the applied uniform field on the scaling behavior 
of the induced gap. 
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